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ABSTRACT 

We investigate the ISM properties of 13 star-forming galaxies within the z ~ 2 COSMOS cluster. We 
show that the cluster members have [N II]/Ha and [O III]/H/3 emission-line ratios similar to z ^ 2 held 
galaxies, yet systematically different emission-line ratios (by ~ 0.17 dex) from the majority of local 
star-forming galaxies. We hnd no statistically signihcant difference in the [N II]/Ha and [O III]/H/3 
line ratios or ISM pressures among the z ~ 2 cluster galaxies and held galaxies at the same redshift. 
We show that our cluster galaxies have signihcantly larger ionization parameters (by up to an order 
of magnitude) than local star-forming galaxies. We hypothesize that these high ionization parameters 
may be associated with large specihc star formation rates (i.e. a large star formation rate per unit 
stellar mass). If this hypothesis is correct, then this relationship would have important implications 
for the geometry and/or the mass of stars contained within individual star clusters as a function of 
redshift. 
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Subject headings: galaxiesistarburst—galaxies:abundances—galaxies:fundamental parameters 


1. INTRODUCTION 

Galaxy clusters formed by the collapse of the highest 
density peak s in primordial density fluctuations (e.g., 
lPeeblesl[l993t IPeacodd 1 199911 . The cluster environment 
strongly affects the star-formation history of cluster 
members. Galaxies in the centers of clusters are ob¬ 
served to be redder and more evolved than galaxies 
on the outskirts of the cluster, as well as field galaxies 


TMl IRaloffh et all 11 9971 1199911 
sive clusters today likely for 
of vigorous star formation ths 
z ^ 1 (IStanford et al.1 11998 

. Thus, the most mas- 
med in short episodes 
it effectively ceased by 
; iBlakeslee et al.l 120061: 

van Dokkum & van der Marell 

20071: Eisenhardt et al. 

2nnSHMei et 31.112009 

:IMancone et al.ll20inl: ITanaka et al. 


cluster formati o n and evolution supp o rt th is scenario 
(iSoringel et al.l 120051 : iDe Lucia et al.l 1200611 . How¬ 
ever, using old stellar populations to reconstruct the 
star formation conditions and to understand why 
star formation ceased in clusters is extremely dif- 
hcult. Recently, several galaxy clusters have been 
discovered in the “redshift desert ” (1.5 < z < 2.5) 
(iKurk et al.l 120091: iPanovich et al.l 120101: iTanaka et alJ 

20101: Fassbend er et al.~ [2Q1I ; iGobat et al. _ 2011: 

Nastasi e t al.1 I20TT: Santos et al.l 20111: iSnitler et ahl 


20121 Stanford et al.l I2012t IZeimann et al.l 120121 : 


Muzzin et al.l 120131 : iT^anaka et al.l 1201311 . With near- 


infrared multi-object spectroscopy on 8-lOm telescopes, 
it is now possible to study the star formation conditions 
in these clusters, prior to (or during) the end of their 
major star-forming episodes. 

Near infrared spectroscopy probes the optical rest- 
frame emission-lines at z 2, providing a wealth of in¬ 
formation about the ionized gas surrounding the young, 
hot stars. The strong lines [N II], Ha, [O III], H/3, and 
the [O II] and [S II] doublets diagnose the power source of 
galaxies, the star formation rate, the gas-phase chemical 
abundance, the ionization state of the gas, the electron 
density of the gas, and the amount o f dust extinction , 
among other physical properties (see lOsterbro"^ I1989L 
for an review). Knowledge of these properties in clusters 
across cosmic time helps to build a more comprehensive 
picture of how star formation progressed in clusters, and 
how the cluster environment impacted the star formation 
process. 

So far, large near infrared spectroscopic studies have 
predominantly focused on field galaxies. High red- 
shift field galaxies appear to have larger [N II] /Ha 

0.3 dex) 


compared with local held galaxies (iHainline et al.l 

2009 

Bian et al. 

120101: IRigby et al.l 12011: Yabe et al.1 

2012 

Erb et al.l 

20061: IZahid et al. 20131: Sha 

□ley et al.1 

2014 

Holden et alJ 120141 ISteidel et al. 12014 

). These ele- 


vated line ratios may indicate a major change in the 
ionized gas or star formation properties as a func¬ 
tion of time. Many possible causes for this off¬ 
set have been proposed, including a larger ionization 
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parameter (e.g., iBrinchmann et all 120081: iShirazi et al.1 
l2013t iKewlev et ^ l2013bll. a harder io nizing radiation 
held (ISteidel et al.1 201^ Stanway et alJ |201-^ . a higher 
ISM pressure (density ) (Shirazi et al. 1^1311. a contri¬ 
bution from an AGN (l&oves et al.l 120061: IWright et al.l 
120101 : iTrump et al.11201 111 and contamination by s hocks 
from galactic outho ws ((Forster Schreiber et al.l 120131 : 
iNewman et al.l (201411 . Sample selection is likely to play 
a large role in the relative importance of each of these 
factors ((Juneau et al.ll201^ . 

If the cluster environment causes more rapid evolution 
during the star-formation stage than in the held, then we 
might expect a difference in the ISM conditions in clus¬ 
ters at high redshift compared with held galaxies. The 
star-formation properties of high redshift cluster galax¬ 
ies relative to the held are currently un der debate. Us¬ 
ing st acked near-infrared spectroscopy, IValentino et al.1 
((2014f) suggested that the relationship between star for¬ 
mation and stellar mass is simi lar among c luster and 
held galaxies a t z ^ 2. However, iTran et ^ (1201011 and 
iBrodwin et al.l ((201311 argue that the star formation rate 
in cluster galaxies is elevated relative to held galaxies at 
the same stellar mass. 

In this paper, we test the hypothesis that star-forming 
conditions in clusters are more ev olved than the hel d 
using the z ^ 2 COSMOS cluster (iSoitler et al.ll20lll . 
The COSMOS cluster is an ideal laboratory for study¬ 
ing the formation of massive nearby galaxy clusters. 
The COSMOS cluster contains four distinct overdensi¬ 
ties ((Spitler et al.ll201^ . Cosmological simulations sug¬ 
gest that at least two of the overdensities will merge 
across cosmic time , creating a Virgo-like cluster by z ^ 0 
((Yuan et al.1120141 1. The cluster has 57 spectroscopi¬ 
cally conhrmed members, and follows the same mass- 
metallicity relation as held galaxies at z 2 (Kacprzak 
et al. 2014, in prep). 

The COSMOS cluster and comparison samples are de¬ 
scribed in Section [2j In Section [3l we describe our near- 
infrared spectroscopy. In Section [4l we compare the opti¬ 
cal line ratios for 10 galaxies in the COSMOS cluster with 
held galaxies at the same redshift, and with local star¬ 
forming galaxies. We investigate the electron density (or 
ISM pressure) in the COSMOS cluster in Section [SJ and 
the ionization state of the gas in Section (6] We dis¬ 
cuss the implications of this work and our conclusions in 
Section [7l Throughout this paper, we adopt the hat A- 
dominated cosmology as measure d by the 7 year WMAP 
experiment {h = 0.72, 0^ = 0.29: lKomatsu et ^1201111 . 

2. SAMPLE SELECTION 
2.1. The z ^ 2 COSMOS Cluster 


The COSMOS cluster was discovered as part 
of the Magellan/FourStar Galaxy Evolution Survey 
(ZFOURG E) photometric red shift survey of the COS¬ 
MOS held ((Scoville et al.l[20?)7l l. The ZFOURGE survey 
uses hve medium bandwidth hlters {Ji,J 2 ,J 3 ,Hs,Hi) and 
the broadband Kg hlter on the Magellan FOURSTAR 
instrument (Perssonl3, Straatman et al. in prep). High- 
redshift galaxy overdensities were identihed using surface 
density maps within narrow Az ^0.2 redshift ranges be- 
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tween 1.5 < z < 3.5. With this technique, the COSMOS 
cluster was revealed with photometric redshifts z ~ 2.2 
(iSpitler et al.ll201^ . Subsequent Keck observations re¬ 
vealed 57 spectroscopically confirmed cluster members 
within a projecte d radius of ^ 5Mpc at z = 2.095 
(jYuan et al.ll2014fl . We compare the ISM properties of 
the COSMOS cluster with two comparison field samples. 

2.2. Local Comparison Field Sample 

We use the Sloan Dig ital Sky Survey (SPSS) as our lo - 
cal comparison sample (|Adelman-McCarthv et ani2006ll . 
Although the COSMOS cluster galaxies will have evolved 
into red elliptical galaxies by z = 0, the SDSS represents 
the range of properties seen in massive (M > IO^Mq) 
nearby star-forming galaxies. W e select galaxi e s from 
the SDSS DR4, as described in IKewlev et al.l (1200611 . 
To avoid aperture effects and incompleteness, we re¬ 
strict the SDSS redshift range to 0.04 < z < 0.1. The 
lower redshift limit corresponds to a minimum aper¬ 
ture co vering fraction of > 20% to minimize aperture 
effects (|Kewlev et al.l 1200a . We use the emission-line 
fluxes from the JHU/MPA catalog with H/3, [O III], Ha, 
[N II], and [S II] emission-lines with signal-to-noise ra¬ 
tios (S/N) > 3tT. These spectra have been corrected by 
the JHU/MPA team for underlying stellar absorption us¬ 
ing stellar population synthesis models. Our SDSS sam- 
ple contains 85,224 galaxies. We remove AGN using the 
IKewlev et al.l (1200611 optical classification scheme, leav¬ 
ing 66,590 star-forming galaxies. 

2.3. z ~ 2 Comparison Field Sample 

We compare the ISM properties of the COSMOS clus¬ 
ter with the properties of 10 serendipitous field galaxies 
between 1.8 < z < 2.2, observed in the same slitmasks 
as the COSMOS cluster candidates. These galaxies were 
identified in the ZFOURGE five medium bandwidth fil¬ 
ters. 

We supplement our small field sample with 39 star¬ 
forming galaxies between 1.8 < z < 2.3 from the Keck 
Baryonic Structure Survey (KBSS-MOSFIRE) survey by 
iSteidel et al.l (|2014ll . The KBSS-MOSFIRE star-forming 
galaxy sample contains ^ 70.9% of galaxies previously se¬ 
lected using rest-frame far-UV spectra, with the remain¬ 
ing galaxies covering a broader range of rest-frame UV 
and optical colors. The final KBSS sample with detected 
[N II]/Ha and [O III]/H/3 emission-lines is not biased in 
star formation rate and stellar mass, relati ve to the par¬ 
ent ( photometric) sample (see Figure 4 in ISteidel et al.l 

[Ml . 

To avoid selection effects based on galaxy mass, while 
ensuring statistically significant samples for comparison, 
we create two comparison field samples, depending on 
the suite of emission-lines available. We refer to these 
two samples as our BPT and [S II] samples, respectively. 

BPT Field Comparison Sample: To compare the po¬ 
sition of galaxies on the standard optical diagnostic dia¬ 
grams (iB aldwin et al.lIlQSii : IVeilleux fc Osterbrocklll98l 
IKewlev et al.l2001atl . we select galaxies with the full suite 

of H/3 A4861, [O III] A5007, Ha A6563, and [N II] A6585 
lines. Steidel et al. apply a lower S/N limit of 5cr to 
the [O III] and Ha emission-lines, and a 2cr limit to the 
[N II] and Ha emission-lines. For consistency, we ap¬ 
ply the same S/N limits to our field galaxies, but we 


highlight galaxies with [N II] or Ha emission-lines with 
S/N< 3cr in our Figures to indicate that these detections 
are marginal. 

The position of a galaxy on the optical diagnostic dia¬ 
grams is sensitive to its gas-phase metallicity, and hence 
to its stellar mass thr o ugh the mass-meta llicitY relation 
(iTremonti et al.l l2004 IZahid et al.l 1201211 . Selecting a 
high redshift sample for the full suite of H/3, [O III], 
Ha, [N II] emission-lines biases a sample towards galaxies 
with detectable [O III] and [N II] lines. This selection can 
introduce a bias against high mass and low mass galax¬ 
ies. Metal-rich nebulae produce low [O III] emission-line 
fluxes and the oxygen lines are sensitive to the electron 
temperature of the gas, which is cool in metal-rich neb¬ 
ulae. On the other hand, selecting for [N II] detections 
biases against low metallicity galaxies which have weak 
[N II] . The [N II] line is sensitive to met allicity only in the 
metel-rich regime (log(0/H)-|-12> 8.4: lKewley fc Dooital 
1200211 because nitrogen changes from a primary to a sec¬ 
ondary nucleosynthetic element at intermediate metal- 
licities. Therefore, we match our field and comparison 
samples in stellar mass. Here, we use a subsample of the 
KBSS-MOSFIRE galaxies matched in stellar mass to our 
ZFIRE cluster sample with 9.0 < M*/M 0 < 10.0. This 
stellar mass selection results in 17 galaxies. 

Our ZFIRE serendipitious field sample contains 4 
galaxies with the full suite of emission-lines H/3 A4861, 
[O III] A5007, Ha A6563, and [N II] A6585 detected to 
signal-to-noise ratios > 3a within the same stellar mass 
range of 9.0 < M*/Mq< 10.0. Applying the KBSS- 
MOSFIRE S/N limits yields one additional galaxy. 

[SII] Field Comparison Sample: We use the [S II] dou¬ 
blet for calculating the electron density of the ionized 
gas in the [S II] zone. To compare the electron density 
of the cluster and field galaxies, we select galaxies with 
rest-frame red emission-lines Ha A6563, [N II] A6585, 
[S II] A6717, and [S II] A6731 detected to the 3cr leyel. 
This selection results in high luminosity, high stellar mass 
samples because the [S II] lines are only weakly depen - 
dent on the metallicity (see e.sf.. lKewley fc DoDitall2003l . 
We therefore match our cluster and field galaxies in stel¬ 
lar mass (9.5 < M*/M 0 < 11.0), resulting in 10 field com¬ 
parison galaxies from our MOSFIRE dataset. Our BPT 
and [S II] field samples are not mutually exclusive, and 
they span the same range of SFRs (6 < < 56). 

For our electron density analysis, we are unable to 
compare with the KBSS-MOSFIRE sample because the 
KBSS-MOSFIRE [S II] line fluxes are not publicly avail¬ 
able. 

3. ZFIRE: HIGH-Z MOSFIRE NEAR-INFRARED 
SPECTROSCOPY 

3.1. Observations 

We obtained near-infrared spectroscopy for the GOS- 
MOS cluster as part of our ZFIRE spectroscopic sur¬ 
vey of high-z cluster candidates. The candidates were 
selected based on t he ZFOURGE photometric redshifts 
([Spitler et al.ll20lll derived from deep images obtained 
through the five ZEOURGE near-infrared medium-band 
filters. The median un certainty for the ZFO URGE pho¬ 
tometry is ~ 0.05 dex ([Tomczak et al.ll^014ll . 

Observati ons of cluster candida tes were made on Keck 
MOSFIRE ([McLean et al.l l2012f) on December 24-25, 
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2013 and February 10-13, 2014. Eight masks were used 
in the iF-band (1.93-2.45 fim) to obtain the Ha and [N II] 
emission-lines, as well as two masks in the i7-band (1.46- 
1.81/rm), to detect the H/3 and [O III] emission-lines. To¬ 
tal on-source exposure times were 2 hours for the K-band 
masks, and 5.3 and 3.2 hours for the two H-band masks. 
Seeing varied between 0.4-0.7 arcseconds throughout the 
exposures. We used a slit width of 0.7”, which gives 
a spectral resolution of i? = 3690 in the K-band and 
R — 3620 in the H-band. The MOSFIRE field of view of 
6.1 X6.1 arcminutes enabled 224 objects to be targeted 
in 6 pointings. We observed an AOV type standard star 
in both the wide-slit mode and the narrow-science-slit 
(~ 0.7” slit width) mode before and after our science 
target exposures for telluric and flux calibration. 

3.2. Data Reduction 

The raw MOSFIRE data were reduced using the 
publicly-available data reduction pipeline (DRP) devel¬ 
oped by the instrument teairQ. The DRP provides 
background-subtracted, rectified and wavelength cali¬ 
brated 2D spectra. All spectra were calibrated to wave¬ 
length in a vacuum. The H-band data were calibrated 
using night sky lines. The K-band data were calibrated 
using both night sky lines and a Neon arc lamp to ensure 
sufficient red emission-lines for wavelength calibration. 
The typical residual for the wavelength solution is < lA. 
We used our own custom IDL routines to correct the 2D 
spectra for telluric absorption, and for flux calibration us¬ 
ing AOV standard stars. We have visually inspected the 
2-dimensional spectra to check for contamination by sky 
lines at the positions of the emission-lines, including the 
weak [N II]A6584 line and the [S II]AA6716,6731 doublet. 
We have also checked that our error spectra fully account 
for skyline residuals. The uncertainty for our flux cali¬ 
bration is ^ 8% (Nanayakkara et al. in prep). The 1-D 
spectra and associated Icr error spectra were extracted 
using an aperture that corresponds to the FWHM of the 
spatial profile. For objects that are too faint to fit a 
Gaussian to the spatial profile, we used the FWHM of 
the stellar profile on the same mask as the extraction 
aperture. 

3.3. Emission-Line Measurements 

Gaussian profiles were fit to the emission lines in H 
and K bands separately. For widely separated lines such 
as [O II]A3727, H/3A4861, single Gaussian functions are 
fit with 4 free parameters: the centroid (or the redshift), 
the line width, the line flux, and the continuum. The 
doublet [O III]AA4959,5007 are initially fit as a double 
Gaussian function with 6 free parameters: the centroids 
1 and 2, line widths 1 and 2, fluxes 1 and 2, and the 
continuum. A triple-Gaussian function is simultaneously 
fit to the three adjacent emission lines: [N II]AA6548, 
6583 and Ha. The centroid and velocity width of the 
[N II]AA6548, 6583 lines are constrained by the veloc¬ 
ity width of HaA6563, and the ratio of [N II]A6548 and 
[N II]A6 583 is const r ained to the theoretical value of 1/3 
given in lOsterbrockI (|1989f l. The line profile fitting is con¬ 
ducted using a minimization procedure which uses the 
inverse of the Icr error spectrum of the DRP. This fitting 

^ See http://code.google.eom/p/mosfire/ 


procedure yields the redshift, line flux, line width, con¬ 
tinuum and the associated statistical errors. Our typical 
3cr flux limit is 1.8 x 10“^® ergs/s/cm^. Example spectra 
of our COSMOS clu ster members are given in Figure 1 
of I Yuan et al.1 (12014( 1. 

From our spectra, we obtained spectroscopic redshifts 
for 180 galaxies. Of these, 57 gal axies are confirmed 
cluster members ((Yuan et al.l l2014l l. A total of 8/57 
ZFIRE cluster galaxies have the full suite of H/3 A4861, 
[O III] A5007, Ha A6563, and [N II] A6585 lines detected 
to signal-to-noise ratios > 3 (t, which we refer to as our 
HPT cluster sample. For consistency with the KBSS- 
MOSFIRE field sample, we also show the positions of 
galaxies with H/3 A4861 and [N II] A6585 lines detected to 
signal-to-noise ratios of 2cr, and [O III] A5007, Ha A6563 
lines detected to 5cr. The KBSS-MOSFIRE S/N selec¬ 
tion criteria would add 5 additional galaxies to the cluster 
sample. We have verified that our results are unchanged 
regardless of whether we apply our S/N > 3 criterion, or 
whether we apply the KBSS-MOSFIRE S/N criteria. 

Our ZFIRE BPT cluster sample was selected to sat¬ 
isfy 9.0 < log(M*/MQ)< 10.0 to match the mass range 
spanned by the BPT field control sample. 

Of the galaxies with Ha A6563, and [N II] A6585 
lines detected to signal-to-noise ratios > 3cr, 4 galax¬ 
ies also have both [S II] A6717 and [S II] A6731 de¬ 
tected to the 3cr level, and span the stellar mass range 
9.5 < log(M*/MQ)< 11.0, where the stellar mass range 
is matched to the [S II] comparison sample. We refer 
to these four [S II] detected cluster galaxies as our [S II] 
cluster sample. Of these, only 3 galaxies have all of the 
H/3, [O III], Ha, [N II], and [S II], emission-lines detected 
to the 3(7 level. 

3.4. Stellar Mass and SFR measurements 

We compute stellar masses using the ZFOURGE 
Jt an ( i AT., photometry. We ht the 

iBruzual fc Gharlc^ (I2003I1 stellar popula tion synthesis 
models with F AST ( Kriek et al.l l2009f l , assuming a 
iGhabrieJ (|2003ll initial mass function, a iGalzetti et al.l 
( 200011 attenuation law, and solar m etallicity. Th e same 
models and parameters were used bv iSteidel et al.l (120141 1 
for the KBSS-MOSFIRE sample. Our BPT cluster sam¬ 
ple, the Steidel BPT field sample, and our ZFIRE BPT 
field sample have the same mean stellar masses, within 
the errors (< log(M*/MQ) >= 9.62 ± 0.06, 9.63 ± 0.08, 
and 9.67±0.07, respectively). Although limited by small 
numbers, our [S II] cluster and held comparison sam¬ 
ples also have consistent stellar masses, within the errors 
(10.0 ± 0.2 c.f. 9.6 ± 0.2, respectively). 

Star formation rates (SFRs) are calculate d from the 
extinc tion-corrected Ha emission-line and the lHao et al.l 
(IMl calibration, as described in detail in Tran et al. 
(2015, ApJ, submitted). We correct the Ha em ission-line 
for extinction using the iGardelli et al.1 (Il989f ) extinction 
curve. The resulting SFRs range between 5 < < 

110 for our BPT cluster sample and 5 < ^ 

for our BPT held sample. The 2D Kolmogorov-Smirnov 
test indicates that our cluster and held samples have sta¬ 
tistically consistent (P(KS)< Icr) extinction and star- 
formation rates, given the sample sizes and errors. We 
do not match the star formation rate range for the clus¬ 
ter and held galaxies because the star formation rate is 
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correlated with the ionization parameter, and the ioniza¬ 
tion parameter is one of the quantities that we aim to 
test for differences among field and cluster galaxies. 

4. OPTICAL EMISSION-LINE ANALYSIS 

4.1. The [N II]/Ha versus [0 III]/Hj3 Diagnostic 
Diagram 

The [N II]/Ha v ersus [O III]/H/3 dia gram was orig¬ 
inally proposed by [Baldwin et al.l (|1981ll to distinguish 
between galaxies powered by H II regions, planetary neb¬ 
ulae, and objects powered by a harder ionizing radiation 
field, such as galaxies containing active galactic nuclei 
(AGN). This diagram (known as the BPT diagram) is 
useful at intermediate and high redshift because the [N II] 
Ha [O III] and H/? lines are often the only emission¬ 
lines observable in near-infrared spectra of high redshift 
galaxies, and the wavelength of the emission-lines are suf¬ 
ficiently close that the [N II]/Ha and [O III]/H/3 ratios do 
not need to be flux calibrated or corrected for extinction. 

Both H II regions and star-forming galaxies 
form a curved locus on the [N II ] /Ha versus 
[O I l l1/H;8 optical diagn ostic diagram (jPoDita et al.l 
120001 : iKewlev et al.l 120061 1. In Figure UKa), we show 
the locus of the star-forming gala, xies in the SDSS sam¬ 
ple. iVeilleux fc OsterbrockI (119^ showed that theoret¬ 
ical photoionization models can reproduce this curved 
locus. This locus is now known as the star-forming abun¬ 
dance sequence because the shape of the locus is dictated 
primarily by the chemical abundance (metallicity) of the 
star-forming regions within galaxies. Metal-poor galax¬ 
ies have low [N II]/Ha and high [O III]/H/3 ratios due 
to the combination of low chemical abundance and the 
electron temperature sensitivity of the oxygen collision- 
ally excited lines. 

The ISM and star formation properties of a galaxy 
affects its position relative to the star-forming abun- 
dance sequence on the BPT diagram (see Figure 1 in 
IKewlev et al.ll 201 ^ . A rise in the ISM pressure causes 
both the [N II]/Ha and the [O III]/H/I line ratios to 
rise. A hard ionizing radiation field (such as from a Wolf- 
Rayet dominated stellar population) causes a similar rise 
in the [N II]/Ha and [O III]/H/3 line ratios. 

However, a rise in the ionization parameter causes the 
[N II]/Ha (and [S II]/Ha) to fall, while [O III]/H/3 rises. 
The ionization parameter is defined as the number of 
hydrogen ionizing photons passing through a unit area 
per second divided by the number density of hydrogen 
atoms uh, and is a measure of the amount of ionization 
that a radiation field is able to produce in an H II region. 
For a spherical geometry, the ionization para meter q can 
be d efined using the Stromgren radius Rg (jStromgrenl 

[HU): 


Qh° 

AttRs’^uh 


( 1 ) 


where Quo is the flux of ionizing photons above the 
Lyman limit. In an ionized nebula, the number density 
of hydrogen is approximately the electron density Ue- 
The ionization parameter is affected by both the hard¬ 
ness of the ionizing radiation field, as well as the bolomet- 
ric luminosity of the ionizing source. The combination 
of theoretical photoionization models, and observations 
of the [O III], II/3, [N II], Ha lines, and the [S II] doublet 


allows one to discriminate between the effect of a dif¬ 
fering ionization parameter, the ISM pressure, and the 
radiation hardness. The [S II] A6717/[S II] A6731 ratio 
constrains the ISM pressure, thanks to the sensitivity of 
this doublet to the hydrogen density of the gas. The 
[S II]/Ha ratio is particularly sensitive to the hardness 
of the ionizing radiation field because the [S II] lines are 
produced in a partially ionized zone that is large when 
the radiation field contains a significant fraction of high 
energy (X-ray and EUV) photons. Finally, the combi¬ 
nation of the [N II]/[S II] and [O III]/H,5 ratios separate 
th e ionization parame ter and the metallicity, as shown 
in iDopita et al.l (|2013ll . 

Figure [IK a) gives the [N II]/Ha versus [O III]/H,5 di¬ 
agram for the SDSS (grey contours) and the COSMOS 
cluster (red symbo ls). The SDSS star-fo rming sequence 
is parameterized bv iKewlev et al.l (|2013afl (their equation 
5). Using this parameterization, we calculate the median 
[O III]/H/3 offset between the COSMOS cluster and the 
SDSS star-forming sequence to be 0.17 dex. We in¬ 
vestigate possible causes for this offset in the following 
sections. 

We compare the position of our COSMOS cluster with 
our field comparison sample in Figure lllb). We have 
matched the stellar mass range of both samples (i.e. 
9.0 < log(M)/M 0 < 10.0) to avoid mass biases. 

Our COSMOS cluster has line ratios that are well 
matched to the field galaxies at the same redshift. 
The two-sided Kolmogorov-Smirnov statistic for the 
[N II]/Ha and [O III]/H/3 ratios yields D-values of 0.34 
and 0.35, respectively, with significance levels of 0.36 and 
0.32, indicating that our cluster and the field sample are 
not significantly different. We have also verified that the 
relationships between [N II]/Ha and stellar mass, and 
[O III]/H,5 and stellar mass are the same for the clus¬ 
ter and field samples, within the errors. We conclude 
that the cluster environment has not produced more 
rapid evolution in the ISM conditions in our star-forming 
galaxies at 2 : ~ 2 . 

4.2. AGN and shock contamination 

AGN strongly affect the [N II]/Ha and [O III]/H/3 op¬ 
tical line ratios. Galaxies containing AGN form a sep¬ 
arate sequence on the BPT diagram, beginning at the 
high metallicity, high mass end of the star-forming abun¬ 
dance sequence (Figured^). This AGN sequence occurs 
because the hard radiation field from an AGN produces 
more ionizations into the O'*''*" ion, and more collisional 
excitations of the [O III] line. The position of an AGN 
host galaxy along the AGN sequence depends on the con¬ 
tribution from star-formation to the line ratios, the slope 
of the EUV radiation field, and the ionizat ion parameter 
of the radiation field ([Groves et al.ll20n4bll . 

We use the redsh ift-dependent theoret ical optical clas- 
sihcation scheme of iKewlev et al.l (l2013bll to identify pos¬ 
sible AGN in our sample (Figured):). This scheme was 
based on the predictions of cosmological hydrodynamic 
simulations, as well as stellar evolutionary synthesis and 
photoionization models. We find that three galaxies in 
our cluster with 2 <S/N([N II] )< 3 may have a small 
(^ 20%) contribution from an AGN based on their opti¬ 
cal line ratios (i.e. [O III]/H/3> 1.0), but these galaxies 
are within the ± 0.1 dex error range of the classification 
line (dotted lines in FigurejS):), and may be low metallic- 
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Fig. 1.— The [N Il]/Ha versus [O IIl]/H/3 diagnostic diagram 
showing our COSMOS cluster (red circles) in comparison to the lo¬ 
cal SDSS sample ( upper panel), our field comparison sample (mid¬ 
dle panel), and the lKewlev et aP (I20131j) redshift-dependent classi¬ 
fication line and error ranges (red solid line and dotted lines respec¬ 
tively, lower panel). The blue lines in the middle panel show the 
position of the local SDSS envelope for comparison. Field galax¬ 
ies are colored according to sample, wh ere pale blue corresp onds 
to the KBSS-MOSFIRE field sample of ISteidel et al.l 1120141) . and 
black corresponds to field galaxies from our MOSFIRE observa¬ 
tions with all strong-lines detected to > 3cr. Orange circles and 
dark blue circles show the additional galaxies that would be added 
to our cluster sample and field sample if we were to apply the S/N 
limits of Steidel et al. 

ity star-forming galaxies with a ha rd ionizing radiation 
field, such as from Wolf-Rayet stars (IKewlev et aUl^Olat 
iDopita et al.ll2006cJ : [Masters et al.ll2014D . or their [N II] 
line fluxes may simply have too low S/N to reliably clas¬ 
sify. We discuss these galaxies further in Section [Si The 
remaining galaxies in our cluster are likely to be purely 
star-forming. 

Shocks can also produce strong [O III]/H^, and 
can mimic AGN on the BPT diagra m locally (e.g., 
IKewlev et al.l I2001bl : IRich et ^ 120111) . At the low 
metallicities of high redshift galaxies, shocks can also 
mimic the emission-l ine ratios of star-forming galaxies 
(|Kewlev et ahll^lSal) . In Figure [2l we show the po¬ 
sition of our COSMOS cluster galaxies relative to the 


slow shock models of IRich et al.l (120101 . 120111) (cyan) 
and the fast shock models of lAllen et al. (I2008D ('green). 
These models are described in detail in IKewlev et al.l 
(|2013al) . and examples of typical mixing sequences are 
given in Rich et al. (2010, 2011). The shock models 
occupy different regions of the diagnostic diagram de¬ 
pending on the metallicity of the sample. Only the low¬ 
est metallicity (log(0/H) -|- 12< 8.17) shock models can 
account for the line ratios of our cluster galaxies. How¬ 
ever, our cluster galaxies lie alon g the standard galax y 
mass-metallicity relation at z ^ 2 (|Kacprzak et al.l2015D , 
with metallicities between 8.3 <log(0/H) -|- 12< 8.7 in 
the IKewlev fc Dooit^ (|2002l) metallicity scale. Here, we 
have c onverted the metalli cities from thelPettini fc Pagell 
(|2004l) calibration into the lKewlev fc Dopital (I2002D scale 
using the calibrations provided in IKewlev fc Ellisonl 
(|2008D to yield consistent metallicity estimates with our 
photoionization and shock models. Because the shock 
models for 8.3 <log(0/H) -|- 12< 8.7 produce optical 
line ratios that are significantly larger than observed (by 
0.1-0.4 dex), we rule out a dominant contribution from 
shocks to the optical emission-line ratios of our cluster 
galaxies. 

Galaxies containing a mixture of ionizing sources, such 
as shocks and star formation will he along mixing se¬ 
quences between the star-forming sequence and the 100% 
shock models. The position and shape of the mixing se¬ 
quences depend on the metallicity of the galaxy and the 
shock velocity. We cannot rule out a small (~10%) con¬ 
tribution from shocks to the emission-lines ratios. High 
angular resolution and high spectral resolution integral 
field spectroscopy with adaptive optics can distinguish 
shock ed regions from s tar-forming regions at high red- 
shift (|Yuan et al.ll2012l) . This integral field spectroscopy 
is one of our future research directions for this sample. 

4.3. Comparison with theoretical models of ISM 
conditions 

We compare the location of our ZFIRE cluster on 
the BPT diagram with our predictions from theoreti¬ 
cal stellar population synthesis and photoio nization mod¬ 
els. Ou r mo dels are descri b ed in detail in IKewlev et al.l 
(l2013al) and iDopita et al.l (120131) . We use t he ch em- 
ical evolution predictions from iDave et'all (|2011l) to 
constrain the gas-phase metallicity history for star¬ 
forming galaxies with M* > 10® Mq . The metal¬ 
licity at a given redshift determines the shape of the 
EUV spectrum produced in the Star burst99 or Pe- 
gase2 evolutionary synthesis models (|Leitherer et al.l 
119991 : iFioc fc Rocca-Volmerangel I1999D . which in turn, 
affects the line intensit ies predicted by the Mappings IV 
photoionization code (IDopita et al.l 1201^ . We use the 
instantaneous and continuous burst models from Star- 
burst99 and Pegase2 with a Salpeter Initial Mass Func¬ 
tion, and the Pauldrach/Hillier stellar atmosphere mod¬ 
els. We note that the choice of IMF makes negligible 
differenc e in the line ratios that we are using in this 
analysis (|Dopita et al.l [201^ . The resulting stellar pop¬ 
ulation spectra are embedded within a spherical neb¬ 
ula. Our Mappings IV photoionization code calculates 
radiative transfer self-consistently throughout the neb¬ 
ula. Mappings IV uses a Kappa temperature distribution 
which is more suitable fo r a turbulent ISM th an a Stefan- 
Boltzmann distribution (|Nicholls et al.l[^12D . For AGN, 
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iGroves et alJ (l2004aD . 

With these models, we produce theoretical predictions 
for how the optical emission-line ratios will appear for 
galaxy samples at different redshifts, based on the fol¬ 
lowing limiting assumptions for star-forming galaxies and 
AGN: 

• Star forming galaxies at high redshift may have 
ISM conditions and/or an ionizing radiation field 
that are either (a) the same as local galaxies {nor¬ 
mal ISM conditions) or (b) more extreme than local 
galaxies {extreme ISM conditions). Extreme condi¬ 
tions in star-forming galaxies can be produced by 
a larger ionization parameter and a denser inter¬ 
stellar medium, and/or an ionizing radiation field 
that contains a larger fraction of photons able to 
ionize 0+ into 0+"'' (i.e. energy > 35.12 eV) rela¬ 
tive to the number of Hydrogen ion izing photons 
(i.e. a harder radiation field). In iKewlev et al.l 
Unni), we use an ionizing radiation field from 
the Pegase 2 stellar population sy nthesis mod¬ 
els (|Fioc fc Rocca-Volmerangjll999ti to mimic the 
hard ionizing radiation field from the stellar at¬ 
mospheres from massive stars with the effects of 
stellar rotation. The difference between our hard 
ionizing radiation field and our normal ionizing ra- 
diation field varies w ith wavelength (see Figure 2 in 
IKewlev et al.ll200Tall . between 1 — lOergs/s/MQin 
log(AFA) between 100 — lOOOA. Similar effects on 
the optical line ratios can be achieved by increasing 
the ionization parameter by up to an order of mag¬ 
nitude (i.e. q = 2 X lO’’ cm/s to g = 2 x 10® cm/s, 
and/or by raising the electron density by up to 
two orders of magnitude (i.e. Ue = 10 cm“® to 
He = 1000 cm“®) (Kewley et al. in prep). 

• The AGN narrow line region at high redshift may 
either (c) have already reached the level of enrich¬ 
ment observed in local galaxies {metal-rich), or (d) 
have the same metallicity as the surrounding star¬ 
forming gas. In the latter case, the AGN narrow¬ 
line region at high redshift would be more metal- 
poor than local AGN narrow-line regions. Our 
AGN narrow-line region models cover a large range 
of power-law indices (—1.2 < a < —2.0) and ioniza¬ 
tion parameters (0.0 < log(17) < —4.0). Therefore, 
the theoretical AGN region may be larger than the 
observed location of AGN on the BPT diagram at 
high redshift. 

These two sets of limiting assumptions yield four limit¬ 
ing scenarios for where galaxies might be located on the 
BPT diagram at different redshifts: 


Fig. 2.— Allen et al. (2008) fast shock models (green) and the 
Rich et al. (2010) slow shock models (cyan) as a function of metal¬ 
licity (panels a-e). The shock model data indicate where galaxies 
that are 100 % dominated by sho cks are likely to lie. The SDSS 
sample from IKewlev et al.l II2006I1 is shown in grey. The positions 
of COSMOS cluster galaxies (red circles) and COSMOS field galax¬ 
ies (black circles) are not consistent with the 100% shock models for 
the metallicities spanned by these galaxies (log(0/H) -|- 12> 8.3). 
Orange and dark blue circles show the additional galaxies that 
would be added to our cluster sample and field sample if we were 
to apply the S/N limits of Steidel et al. 

we use the dusty radiation-pressure dominated models of 


1. Normal ISM conditions, and metal-rich AGN 
narrow-line regions at high-z. 

2. Normal ISM conditions, and metal-poor AGN 
narrow-line regions at high-z. 

3. Extreme ISM conditions, and metal-rich AGN 
narrow-line regions at high-z. 

4. Extreme ISM conditions, and metal-poor AGN 
narrow-line regions at high-z. 






















Our theoretical model predictions for each of these four 
limiting scenarios are given by the solid lines in Figure |3] 
for z = 0 (upper panel) and z 2 (lower panel). The 
COSMOS cluster appears to span a range of ISM and 
star-forming conditions, from nearly local-star-forming 
conditions (panels 1, 2) to more extreme star-forming 
conditions (panels 3, 4). Our data do not contain suffi¬ 
cient AGN to distinguish between the two limiting AGN 
scenarios. We conclude that at least 4 of our cluster 
galaxies have more extreme star forming conditions than 
local star-forming galaxies, on average, and that their 
star-forming conditions are similar to field galaxies at 
the same redshift, as seen in Figure 1. 

5. ELECTRON DENSITY/PRESSURE OF THE ISM 

We investigate whether the [N II]/Ha ratio correlates 
with the electron density of the ionized gas. The mean 
ISM pressure, P, is related to the total gas density, n, 
and mean electron temperature, Tg, through n = 

Here, the total density n is related to the electron density 
through n = 2ne(I + H^,/H). In a fully ionized plasma, 
the electron temperature is ~ 10^ K and the ISM pres¬ 
sure is directly proportio nal to the electron density (see 
e.g., iDooita et al.l 1200^ for a discussion). The mean 
ISM pressure is related to the mean mechanical luminos¬ 
ity flux from the central star clusters within H II regions. 
In the standard adiabatic shell model for H II regions, 
H II regions expand until the internal pressure from 
mass-loss and supernova energy- driven bubbles equal s 
the ambient pressure of the ISM (|Oev fc ClarkdFlQQTl l. 
Therefore, in the expanding bubble model, the ISM pres- 
sure is inversely pro portional to the H II region radius 
(|Dopita et al.ll2006bll . 

We calculate the electron density using the [S II] dou¬ 
blet ratio, [S II] A6717/[S II] A6731. We assume a four- 
level model atom in our Mappings IV models. The rela¬ 
tionship between the [S II] A6717/[S II] A6731 and elec¬ 
tron density given by Mappings IV is identical to th e 
[S II]-electron density relationship in lOsterbrockI (II989D . 
When the [S II] ratio approaches ~ 1.45, the electron 
density is in the low density limit (i.e., < 10 cm~^). 
Only one out of five of the COSMOS cluster galaxies 
with measurable [S II] ratios is in the low density limit. 

In Figure S] we show the electron density versus the 
[N II]/Ha ratio for the COSMOS cluster and the z ^ 
2 field galaxies from our [S II] sample. The cluster 
galaxies span a similar range of electron densities as 
the field galaxies at the same redshift. The two-sided 
Kolmogorov-Smirnov statistic yields a D-value of 0.4, 
with a significance level of 0.54, indicating that the two 
samples are not significantly different. Clearly, there is 
no systematic difference in the electron density among 
our cluster and field galaxies at z ~ 2, within the limita¬ 
tions of our small sample size and relatively large errors. 

6. THE IONIZATION STATE OF THE GAS 

The ionization state of the gas (or ’’ionization pa¬ 
rameter”) is traditionally measured using the ratio of 
two ionization states of the same atomic species. The 
most commonly used ionization parameter diagnostic 
is the [O III] A5007/[0 II] A3727 ratio. This ratio 
is usually calibrated as a function of ionization pa¬ 
rameter and metallicity using stellar population synthe- 
sis and photoionization models (e.g., iKewlev fc Dopital 


l2002t iKobulnickv fc KewlevI 1200411 . Here, we use the 
[N II]/[S II] and [O III]/H/I ratios to distinguish the 
ionization parameter from metallicity, because the rest- 
frame blue [O II] line is not available within the wave¬ 
length coverage of our spectra. The [N II]/[S II] ratio is 
sensitive to ionization parameter because the ionization 
potentials of the [N II] A6584 and [S II] AA6717,31 lines 
differ significantly. The ratio of [N II]/[S II] is also sen¬ 
sitive to metallicity, particularly for metallicities below 
log(0/H) -I- 12< 8.69. This metallicity sensitivity occurs 
because the ratio of collisional excitation rates of [N II] 
and [S II] is a weak function of nebular temperature, 
which becomes high at low metallicities due to the lack 
of coolants in the nebula. We note that although sulphur 
is a typical alpha element (i.e. its abundance sc ales di¬ 
rectly proportional to oxygen) (|Rvde et al.ll2014fl . nitro¬ 
gen has both primary and secondary nucleosynthetic ori¬ 
gins, and is sensitive to the metallicity in the secondary 
nucleosynthetic regime (log(0/H) -|- 12> 8.4). To ac¬ 
count for this effect, our models use an empirical function 
for N/0 versus 0/H which takes into account both the 
primary and secondary nucleosy nthetic components, as 
described in IDooita et al.l (|2013[ 1. The [O III]/H/3 ratio 
helps to separate the ionization parameter from metal¬ 
licity, particularly below log(0/H)-|-12< 8.69, where the 
[N II]/[S II] ratio becomes insensitive to ionization pa¬ 
rameter. In this regime, the [O III]/H/3 ratio is primarily 
sensitive to the ionization parameter. At higher metallic¬ 
ities, log(0/H)-|-12> 8.69, the oxygen lines in the optical 
and the infrared dominate the cooling of the nebula, lead¬ 
ing to less collisional excitations of the 0+''' ion, and a 
stronger dependency between [O III]/H/3 and metallicity. 

In Figure [S] we show theoretically how the ioniza¬ 
tion parameter (yellow-red curves) and the metallic¬ 
ity (green-blue curves) are related on the [N II]/[S II] 
versus [O III]/H/3 diagram. For ionization parameters 
log(q) < 7.75, the [N II]/[S II] and [O III]/II/3 line ra¬ 
tios successfully separate the ionization parameter from 
metallicity. We show cluster and z ~ 2 field galaxies 
with > 3cr detections of [SH] (red circles and black cir¬ 
cles respectively), as well as cluster galaxies with less 
reliable 2cr detections of [SH] (orange circles). The ma¬ 
jority of local star-forming galaxies have ionization pa¬ 
rameters between 6.0 < log(g) < 7.25, while our COS¬ 
MOS and field galaxies have larger ionization param¬ 
eters, i.e. log((?) > 7.25. The two galaxies with the 
largest [O III]/H,5 ratios lie above the theoretical curves, 
suggesting a possible small (i.e. ~ 20%) AGN contri¬ 
bution. Additional information from the X-rays, IR, or 
radio is required to confirm the presence of an AGN in 
these galaxies. 

The electron de nsity of our models (described in 
IDooita et al.l l20I4jl is fixed by the ISM pressure, and 
corresponds to Ue ~ 10 cm“^. An electron density of 
10 cm“^ is consistent with all of our cluster galaxies and 
2/3 of our field galaxies shown on Figure [H within the 
errors. We note that if the model electron density is in¬ 
creased from 10 cm“^ to an extreme value of 1000 cm“^, 
the theoretical [O III]/H/3 ratio would rise by ^ 0.4 dex. 
The [S II] line fluxes are subject to large uncertainties 
(up to 1000 cm“^). Within these large errors, two clus¬ 
ter galaxies and one field galaxy could potentially have 
electron densities as large as 1000 cm“^. We highlight 
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Fig. 3. — The [N Il]/Ha versus [O IIl]/H/3 diagnostic diagram showing the local SDSS galaxies (top panel), and our z = 2.095 COSMOS 
cluster (bottom panel; red circles). Orange and dark blue circles show the additional galaxies that would be added to our cluster sample and 
field sample if we were to apply the S/N limits of Steidel et al. Solid lines show our theoretical predictions for the star-forming abundance 
sequence (left curves) and the starburst-AGN mixing sequence (right curves) for z=0 (top) and z=2 (bottom) for 4 limiting model scenarios 
(columns). Column (l): Normal ISM conditions, and metal-rich AGN NLR at high-z; Column (2): Normal ISM conditions, and metal-poor 
AGN NLR at high-z; Column (3): Extreme ISM conditions, and metal-rich AGN NLR at high-z; Golumn (4): Extreme ISM conditions, 
and metal-poor AGN NLR at high-z. Blue dotted lines in the lower panel show the position of the local models, for comparison. 



Fig. 4. — The [N II]/Hq: ratio versus electron density for the 
COSMOS cluster (red circles), and field galaxies observed as part 
of the ZFIRE survey (black circles). Also shown are additional 
galaxies that would be added to our cluster sample and field sample 
if we were to apply the S/N limits of Steidel et al. (orange and 
dark blue circles, respectively). 

the position of these three galaxies on Figure [S] with a 
circular outline. Even with an 1000 cm“^, two out of 
the three high electron density galaxies would still have 


ionization parameters larger than local SDSS galaxies at 
the same stellar mass. 

Large ionization parameters have been seen 
in high redshift galaxies in previous work (e.g.. 

Brinchmann et al.1l20n8l : iShirazi et al.ll201.‘il : lSteidel et alJ 

2014fl . but t he cause of such a large ion izati on parameter 
is unk nown. iBrinchmann et all (120081 1 and IShirazi et al.) 
(|2013ll suggest that higher electron densities and a larger 
escape fraction of H-ionizing photons may be responsible 
for t he larger ion i zation parameter. On the other 
hand, iSteidel et al.l (|20l4 l suggest that a high effective 
temperature in the massive stellar population, possibly 
from stellar rotation and/or binaries may be responsible 
for the large ionization parameters. Our sample is too 
small to reach statistically significant conclusions on the 
cause of the high ionization parameter in our cluster and 
field galaxies, especially given that four of our galaxies 
lie along the edge of the ionization parameter-metallicity 
grid in Figure [5] where distinguishing between ionization 
parameters above log( 9 ) > 7.75 is impossible. We 
note that the cluster galaxy with the lowest ionization 
parameter {log{q) ^ 7.4) has a significantly lower 

specific star formation rate (log(SFR/M,) ^ —8.9) 
than the majority of our cluster galaxies (mean 
< log(SFR/M,) >= —8.4 ± 0.2). If there is a relation¬ 
ship between specific star formation rate and ionization 
parameter, then this may have important implications 
on the geometry and mass of the star-forming clusters 
at high redshift. 

7. CONCLUSIONS 

In this paper, we have investigated the rest-frame opti¬ 
cal line ratios of the COSMOS galaxy cluster at redshift 
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Fig. 5.— The [N Il]/[S II] ratio versus [O III]/H^ for the 2 : ^ 2 
COSMOS cluster for > 3 (t detections of [SIl] (red circles) and 
for ~ 2cr detections of [SII] (orange circles), compared with our 
2 : ~ 2 ZFIRE field galaxies (black filled circles) and the local 
SDSS sample (dots). The stellar mass range of the SDSS sam¬ 
ple has been matched to the stellar mass range of our z 2 
galaxies (9.5 < log(M)/M 0 < 11.0). The colored curves show 
our theoretical photoionization models for a range of metallicity 
(blue-green lines), and ionization parameter (yellow-red lines), as 
indicated in the legends. The COSMOS cluster and field galaxies 
have larger ionization parameters than the majority of local SDSS 
galaxies. Galaxies with electron densities that could be as large as 
Tie ~ 1000 cm“^ within their errors are highlighted with a bold 
circular outline. 

2.095. We compare the line ratios of the COSMOS clus¬ 
ter with field galaxies at the same redshift. We incor¬ 
porate theoretical stellar evolution, photoionization, and 
shock models into our analysis to investigate any changes 
in ISM conditions, ISM pressure, or ionization parameter 
as a function of redshift or as a function of environment. 

We find that the COSMOS cluster optical line ratios 
are indistinguishable from field galaxies at the same red¬ 
shift. We find no statistical difference among the ISM 
pressure or ionization parameter of field or cluster galax¬ 
ies, within the limited size of our sample. We conclude 
that for the COSMOS cluster, the cluster environment 
makes no or limited impact on the ISM properties of its 


member galaxies. We also rule out shocks as a dominant 
contributor to the optical emission-lines in our cluster 
galaxies. 

We show that the COSMOS cluster has larger 
[N II]/Ha and/or [O III]/H/3 ratios than local galaxies, 
similar to the large line ratios seen in field galaxies at 
high redshift. This offset is consistent with a theoretical 
change in ISM conditions with redshift. We investigate 
possible causes of this change of ISM conditions using our 
theoretical models. We rule out a change in ISM pressure 
as the dominant cause of the extreme ISM conditions in 
our high redshift sample. We use the [N II]/[S II] versus 
[O III]/H/3 diagram to distinguish the ionization param¬ 
eter from metallicity. We show that the large line ra¬ 
tios are likely to result from a large ionization parameter 
(log((;) > 7.25). The ionization parameters in our COS¬ 
MOS cluster and field galaxies are significantly larger (by 
up to 1 dex) than seen in local galaxies. We suggest that 
there may be a relationship between ionization param¬ 
eter and specific star formation rate. A larger sample 
of high-redshift galaxies with ionization parameter and 
specific star formation rate measurements is required to 
test this idea and to explore possible implications on the 
H II region geometry and star cluster properties. 
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